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Abstract 

We address the problem of the apparently very small magnitude of CP violation in the stan- 
dard model, measured by the Jarlskog invariant J. In order to make statements about probabili- 
ties for certain values of J, we seek to find a natural measure on the space of Kobayashi-Maskawa 
matrices, the double quotient U(1) 2 \SU(S)/U(1) 2 . We review several possible, geometrically 
motivated choices of the measure, and compute expectation values for powers of J for these mea- 
sures. We find that different choices of the measure generically make the observed magnitude of 
CP violation appear finely tuned. Since the quark masses and the mixing angles are determined 
by the same set of Yukawa couplings, we then do a second calculation in which we take the 
known quark mass hierarchy into account. We construct the simplest measure on the space of 
3x3 Hermitian matrices which reproduces this known hierarchy. Calculating expectation values 
for powers of J in this second approach, we find that values of J close to the observed value are 
now rather likely, and there does not seem to be any fine-tuning. Our results suggest that the 
choice of Kobayashi-Maskawa angles is closely linked to the observed mass hierarchy. We close 
by discussing the corresponding case of neutrinos. 

PACS numbers: 12.15.Hh, 14.60.Pq, 02.20.Hj, 02.40.-k 
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1 Introduction 



A traditional attitude to theoretical physics has been that the main problem is to discover the 
fundamental laws of physics and leave it to experiment and observation to decide what particular 
implementation best describes "Our Universe." Thus traditionally a "physical theory" is often 
thought of in terms of a local Lagrangian including certain "coupling constants," "mass ratios," and 
"mixing angles," all of which, since Planck's introduction of Planck units [1], may be taken to be 
dimcnsionlcss numbers. In addition, the local Lagrangian must be supplemented with an account of 
the general class of boundary conditions for which the variational principle is valid. Different classes 
of boundary conditions are usually thought of as different "superselection sectors" of the theory, and 
describe qualitatively different types of situations which traditionally are not thought of as having 
any relation to one another. 

Within each sector, there are many solutions of the equations of motion, each of which may be 
specified by providing suitable "initial conditions." Classically these conditions may be thought of 
as the space of classical histories, and given in terms of Cauchy data modulo the relation that two 
sets of Cauchy data giving the same history are taken to be equivalent. Quantum mechanically one 
thinks in terms of some initial, and thus in the Hciscnbcrg picture, eternal state. 

The hope has frequently been expressed in the past that eventually theorists will hit upon a 
unique theory, with all coupling constants determined by consistency or symmetry considerations, 
and with just one superselection sector. Even given such a "theory of everything" (TOE), there 
remains the issue of boundary conditions or initial state, as emphasized by Hawking [2]. Recently, 
however, there has been a considerable decline in optimism and few now seem to believe in a single 
TOE with a single superselection sector, and many refer to a "landscape" of theories. 

One approach to this perceived crisis in theoretical physics is to resort to "anthropic" considera- 
tions and invoke the idea that there may indeed exist, in the Platonic sense, an enormous number of 
"possible universes," of which only very few will allow the development of sentient beings, and even 
fewer will allow sentient beings like ourselves. Thus one is led to contemplate the ensemble of all pos- 
sible universes, sometimes referred to as a "multiverse" [3]. This ensemble is sometimes thought of 
non-Platonically^ as an ensemble of connected subsets of a much bigger physically existing universe, 
referred to as a meta- universe [4] . 

At this point it may be helpful to remark, lest the daunting task of thinking about and making 
more precise, the nebulous notion of such a multiverse should not be thought entirely a problem for 
theorists seeking credit for making predictions about the world we see about us: that the observers 
and experimenters must also face up to that task when assessing the reliability of their measurements 
or the extent to which they can confirm theoretical predictions. All such activities are essentially 
Bayesian in character [5] and require some notion of "priors," that is, some sort of a priori measure 
of the space of possibilities. 

This problem has been addressed, with admittedly only partial success, in a previous paper [6] 

1 in an Aristotelian sense? 
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where the multiverse, for concreteness, was identified with the set of classical histories of a min- 
isuperspace cosmological model. A well-defined and natural local measure on the space of classical 
histories is easily constructed, but unfortunately the total measure of all histories, even in this finite- 
dimensional truncation of the full set of solutions of Einstein's equations, is infinite. The problem 
was recently revisited in [7]. 

In the present paper, we shall turn to the problem of finding a natural measure on the space of 
coupling constants. Thus the multiverse in the present paper is a set of Lagrangians parametrized 
by a manifold X or "moduli space," whose coordinates consist of masses, mixing angles, coupling 
constants, etc., and we wish to place a natural measure on this space. We hope this will be useful 
for anthropic considerations such as those of [8], where X = S , the circle parametrizing the phase 
of the axion. In that case the issue of a measure was trivial, but in more complicated cases such 
as we shall consider in the present paper, the situation is more complicated. We also hope that 
the work in this paper will help in clarifying the notion of "fine-tuning," which is so prevalent in 
phenomenological discussions. 

The structure of this paper is as follows. After introducing the notion of geometric probability and 
outlining the Kobayashi-Maskawa theory of CP violation in the standard model, we discuss metrics 
on SU(3) and its quotients in Sec. [2J starting with a left-invariant metric on SU(3) which induces 
a metric on the flag manifold SU(3)/U(l) 2 . We perform a Kaluza-Klein type reduction on the left 
phases and discuss different possible metrics on the double quotient U(1) 2 \SU(3)/U(1) 2 , the space 
of Kobayashi-Maskawa matrices. We also discuss the metric used by Ozsvath and Schiicking [9], 
and argue that it lacks a geometrical justification. 

We then use all metrics we have discussed to compute statistics of the Jarlskog invariant J. While 
the measure on SU(3)/U(1) 2 is independent from the choice of left-invariant metric, the measure 
on the double quotient is nonunique. We find that in each case the standard deviation AJ (with 
(J) = 0) is about three orders of magnitude greater than the experimentally observed value of J, 
which appears to be finely tuned. In Sec. [H we do a closer numerical analysis of the probability 
distribution of | J| on the double quotient, using several possible choices for the measure. We quantify 
the statement that a magnitude of CP violation as small as observed appears unlikely. 

In Sec. [51 we take a different viewpoint: We now consider random distributions on the space of 
mass matrices in the standard model. We therefore need to find a measure on the space of 3 x 3 
Hermitian matrices. We find that the simplest choice which gives convergent integrals over this 
space, and expectation values for squared quark masses which reproduce the observed values, is a 
Gaussian weighting function with four free parameters, which can be chosen appropriately. We then 
find that the standard deviation AJ is much smaller for this measure, making the observed value 
of J appear typical. We conclude that with an appropriate distribution which respects the known 
quark mass hierarchy there is no need for fine-tuning in J. 

We briefly discuss the case of neutrinos in Sec. El explaining the general theory, and the difference 
between Dirac and Majorana masses. We cannot give reliable predictions for A J for neutrinos, due 
to the absence of known values for their masses. 
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1.1 Geometric Probability 

The construction of appropriate measures over spaces of geometric objects goes back to the 18th 
century cosmologist Buffon and his celebrated needle problem [10]. The reader may find a general 
account of the subject in [11]. The simplest case to consider is when the space of coupling constants 
X may be regarded as a finite-dimensional homogeneous space with respect to some Lie group G of 
symmetries, and the stabilizer or little group is H C G. Thus X = G/H. If dimX = n, our aim 
is to construct an n-form on X which is invariant under the action of G. In the case that X = G, 
as in the example of the axion circle, this is completely unambiguous. We pick any n-form at the 
unit element e e G and spread it over G by left or right translation. On a unimodular group, such 
as a compact group or a semisimplc group, left or right translation will give identical results. The 
original n-form, being a top degree form, is unique up to a multiple. This multiple can be fixed by 
normalizing the total measure to unity. The normalized measure is therefore unique. 

We could, if we wished, construct the measure as the Riemannian volume element of any left or 
right-invariant metric on G. The result would be the same. In practice, a convenient procedure for 
calculating the measure could be to construct an invariant metric on X = G and then calculate its 
Riemannian volume element. Often, the bi- invariant or Killing metric is the most convenient choice. 

In the case of a coset, X = G/H, the measure can again be taken to be any n-form at some 
arbitrarily chosen point x e X, which is then spread around using the group action. Since any 
n-form at x will be H- invariant, the result is again unique and invariant under all the symmetries 
of the problem. Of course it is possible that one may express X = G/H in more than one way. 
This could in principle give rise to some discrete nonuniqueness, but in practice this seems not to 
be important. 

Although the situation when coupling constants may be regarded as belonging to a homogeneous 
space is quite satisfactory, it is often the case that coupling constants belong to an inhomogeneous 
space. In particular, in the case of "mixing angles," they typically belong to a double coset, or 
bi-quotient, of the form Hi\G/H 2 , where H\ and H2 are (not necessarily identical) Lie subgroups 
of G. The reason for this is that mixing angles relate two unitary bases for the same space of 
physical states. The two unitary bases may not be unique. In particular, it is often the case that 
the individual basis vectors can be multiplied by arbitrary phases. In this case, Hi and H 2 may 
belong to U(l) k , where k is the number of states in the basis. In the case of the Kobayashi-Maskawa 
matrix the states are quarks, and one basis diagonalizes the strong Hamiltonian while the other 
basis diagonalizes the weak interaction quantum numbers. 

A biquotient, or double coset, Hi\G/H 2 , is typically not a homogeneous space. This is because 
the left action of G will not in general commute with Hi, and similarly, the right action of G 
will not commute with H 2 . As a consequence, one cannot, in the case of biquotients, use group 
invariance to construct an unambiguous measure on the space of mixing angles. Later in this paper, 
we shall explore in detail some available options, and the extent to which they affect the probability 
distribution of mixing angles. 
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1.2 CP Violation and the Jarlskog Invariant 

In this section we shall review the Kobayashi-Maskawa theory of CP violation in the quark sector 
of the standard model. 

If m and m! are the (Hermitian) mass matrices for the charge | and — ^ quarks, respectively, 
then there exist unitary matrices U and U' such that 

UmU^ = diag(m«, m c , mt) , U'm'U — diag(m<j,m s ,m6) . (1) 

The Kobayashi-Maskawa matrix V is defined by 

V = UU' f . (2) 

The normalized mass eigenstates are only defined up to a phase, and changing these phases changes 
the matrices U and V according to 

U^P L U, V—tP^U'. (3) 

Hence the Kobayashi-Maskawa matrix changes according to 

V -> PlVPr , (4) 

where Pl and Pr are diagonal matrices belonging to SU(3). In other words, Pl and Pr may each 
be thought of as belonging to T 2 = U(l) x t/(l), the maximal torus of SU(3). Thus the four- 
dimensional space of CP violating parameters should be thought of as an element of the double 
coset, or biquotient, U(1) 2 \SU(3)/U(1) 2 , whereas the matrices U and U' should be thought of as 
elements of the left coset U(l) 2 \SU(3). 

In the discussion of geometric probability attempted in this paper, one could take the viewpoint 
that U and U' are the fundamental objects relevant in CP violation, which would lead to discussing 
distributions on (U (1) 2 \SU (3)) 2 . One can then use the fact that only V = UU'^ appears in the 
Kobayashi-Maskawa theory to reduce this to a distribution on a single U(1) 2 \SU(3), as we shall see 
in Sec. 13.11 Alternatively, one considers V as fundamental and considers the biquotient. 

Because the right action of U(l) 2 is free, the intermediate coset SU (3)/C/(l) 2 is a compact smooth 
homogeneous space without boundary, on which SU(3) acts by left actions. In fact SU(3)/U(1) 2 
is an example of a flag manifold. The maximal torus U(l) 2 acts on the flag manifold via left 
actions of SU(3), but its action on SU(3)/U(1) 2 is not free, and as a consequence, the biquotient 
C/(l) 2 \5'C/(3)/J7(l) 2 is not a smooth compact manifold without boundary. Rather, it is a stratified 
set whose boundary consists of components at which either or both of the left-acting U(l) factors 
has fixed points. 



In the standard notation 



^V u d V us V u b\ 
V= V cd V cs V cb , (5) 
\V td V ts V tb j 
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and it is customary to choose the phases so that 

(l W C13 

C 23 S 23 
\0 -S23 C 23 / 



si 3 e- 4 A 



\-Sl3 



1 

s 



-i<5 



Cl3 / 



' C12 8X2 (A 
-Sl2 C12 

V 01/ 



(G) 



where s\ 2 — sin#i2, c\2 — cos#i2, etc., and the angles 6\ 2 , #13, and #23 are taken all to lie in the 
first quadrant {i.e. between and ^7r). 

One can take the angle 5 as a measure of CP violation, but its definition depends on the choice 
of phases. Jarlskog [12, 13] introduced a formalism that eliminates this arbitrariness. She defined a 
Hermitian tracefree matrix C by 

[m, m] = iC , (7) 



and took det C as a measure of CP violation. She showed that 



where 



detC = -2TB J , 



[m t - m u )(m t - m c )(m c - m u ) , B = (m b - m d )(m h - m s )(m s - m d ) , 



and the Jarlskog invariant J is given by 



J = 3m(V 11 V 22 V 1 * 2 V 2 * 1 ) 



(8) 



(9) 



(10) 



Despite appearances, J is independent of the arbitrary phases. In other words, it is invariant 
under In fact, it has an extremely elegant geometrical interpretation. Since V is a unitary 
matrix, its three rows and columns are orthogonal. Thus, for example, there are three relations of 
the form 

{VV% 2 = V n V 2 \ + Vx2 V 2 * 2 + V 13 V. 2 \ - . (11) 

The three complex numbers a — Vn V 2 \, b = V\ 2 V 22 , and c = V13 V 23 , satisfying a + b + c = 0, may 
be thought of as the three sides of a unitarity triangle in the complex plane. The absolute value of 
J is twice the area of this triangle: 



|J| = \3m(ab*)\ = \3m(ac*)\ = \3m(bc*)\ . 



(12) 



The effect of the transformation Q is to rotate this triangle in the complex plane, but the area 
\\J\ is unchanged. Less obviously, the same area results from taking either of the two other possible 
inner products, (FT^t) 13 = or (VV^) 2 3 = 0. Thus J is an invariant, and so it is well defined on the 
space of mixing angles. 

In terms of the standard parametrization ©, the Jarlskog invariant is given by 



J = C12 c 23 c 2 13 S12 s 23 si 3 sin <5 . 



(13) 
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One could choose to take a different quantity as a measure of CP violation. Jarlskog [14] suggested 
appropriately normalizing the determinant ([5]) and using 

detC 



acp 



3-y/6- 



(14) 



(TrC 2 ) 3 / 2 ' 

which takes values between +1 and —1 and is zero if and only if CP is conserved. Written out 
explicitly in terms of the quark masses and mixing angles, this is a complicated expression that we 
do not give here. As in the present paper the observed quark mass hierarchy is assumed, we shall not 
consider the case of coinciding quark masses, and we concentrate on J as a measure of CP violation. 

Another possible source of confusion is the assumption of general, not necessarily Hermitian, 
mass matrices. In this case the commutator ([7]) is replaced by 



mm* ,m'm'^ — 



]=iC 



(15) 



in order for C to be Hermitian. The use of C or C can lead to ambiguous "orders of magnitude" 
estimates for CP violating processes, e.g. when discussing baryogenesis. We will assume that m and 
m' are Hermitian, and as our calculations only involve J these considerations will not be relevant. 



2 Metrics on SU(3) and its Quotients 

A generic element U of SU(3) is conveniently parametrized by eight real coordinates (p, q, r, t, x, y, z, w), 
so that 

U = T L WT R , 



where 



and 



with 



Tr = e 



W — e ixXr e~ iwX3 e lyX -' e iwX3 e izX2 



0<X<\lt, 0<J/<±7T, 0<Z<i7T, 0<w<2tt. 

Here, we are using the standard Cell-Mann representation for the generators of SU(3): 



A a = 



A, = 



/o 1 o\ 

10 0, 

V0 Oy 

(° ° ~^ 



V o/ 

A o o\ 

0-10 

^0 Oy 



A 2 = 



A fi = 



fo -i o\ 

i , 

yo o 0y 

/o o o\ 

1, 

V° 1 °/ 

(l \ 

1 

v -2 j 



A 4 = 



X 7 = 



(o 





A 











V 







(o 





°) 








— i 


V° 


i 





1 

71 



(16) 
(17) 
(18) 
(19) 



(20) 
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Explicitly, the matrices W, Tl, and Tr are given by 



W 



CyC Z 



^x J y^-z 
-'X $y Qz 



S T C 



CyS Z 
— 6 St SyS 



X^Z 



x J y J z 
xSyS z 



Oy 

C-xCy J 



T L = disig(e 2ip ,e- ip+iq ,e- ip - iq 



liag^e 



diag( e ir+lt , e ir - i{ , e- 2ir 



(21) 

(22) 
(23) 



where we use the notation s x = sinx, c x — cosx, etc. 

If we identify W as the Kobayashi-Maskawa matrix V in the standard conventions I©, then 



X = 023, V = 0l3, Z = 6 12 , 

If we define left-invariant one-forms a a by 

U- 1 dU = i\ a a a , 
then the general left-invariant metric on SU(3) takes the form 

ds 2 = 9ab Co Ob , 



w 



(24) 



(25) 



(26) 



where g a b is a constant symmetric matrix. 

For a general choice of the matrix g a i,, the metric admits no further isometries beyond the left 
action of SU(3), which we denote by SU(3)l- For special choices of g a b, however, the metric is 
additionally invariant under the right action of some subgroup K of SU(3)r. The most symmetric 
such case, the bi-invariant or Killing metric for which K is the full right-acting SU(3)r, arises if 
g a b is proportional to 8 a b- The various intermediate possibilities, of which there are five, are listed 
in [15]. 

In the generic case (i.e. when K is the identity), 28 = 36 — 8 parameters are required to specify 
the metric. One of these parameters sets the overall scale. For the intermediate cases there are 
correspondingly fewer parameters [15]. The bi-invariant metric has the smallest number, namely 
just the overall scale. In all cases, the invariant measure on the group SU (3) is the same and given 
by 

i x = N\\a a , (27) 



where TV is a constant normalization factor. 

One of the intermediate cases given in [15] corresponds to 



ds 2 



a 2 (a\ +aj+ a 1 ) + b 2 (a\ + a\ + a 2 ) + a 2 (o| + erf) . 



(28) 



This has the symmetry SU(3)l X 5*0(3)^, where the SO(3)r is generated by A2, A5 and A7. Re- 
markably, there is a second Einstein metric in this class [16], in addition to the standard bi-invariant 
metric that arises when a = b. The nonstandard Einstein metric occurs when b = a/\/TT. 
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However, the measure on the biquotient U(1) 2 \SU(3)/U(1) 2 of SU(3) is not unique. In partic- 
ular, if one constructs the measure from an invariant metric it will depend upon the metric that 
is used. One of the cases enumerated in [15], which is of particular interest for purposes, is when 
K = U(l) x U(l). One may denote these U(l) subgroups by £^(1)3 and U(l)s, indicating that they 
are generated by A3 and As- The possible U(l)s x U(l)$ x SU(3)l invariant metrics on SU(3) are 



ds 2 =01(0-1 + a l) + P {°l + °i) + 7 (^6 + 07) + <*l °3 + ^2 <4 + 25 3 cr 3 a s . 
The induced metric on the right coset SU(3)/U(1) 2 is then given by 

ds 2 =a(a 2 + of) + p {p\ + of) + 7 (oe + • 
The normalized invariant measure on this coset is given by 

/i = N <j\ A 02 A 04 A 05 A (76 A (77 . 



(29) 



(30) 



(31) 



There is no similarly unique construction of a measure on the biquotient J7(l) 2 \5E/(3)/f7(l) 2 , 
because there is no natural action of SU(3)l on it. The reason for this is that the U(l) 2 of the left 
quotienting is the maximal torus in SU(3)l, and so no other generators commute with it. 

Locally, the biquotient U(1) 2 \SU(3)/U(1) 2 is a fiber space whose fibers are orbits of U(1) 2 L x 
U(1) R , whose action has fixed points. The biquotient is therefore not a smooth manifold. Never- 
theless, any metric on SU (3) will induce on any local section a metric, and hence a Riemannian 
measure. However, the metric and the measure will in general depend upon the choice of section. In 
the language of Kaluza-Klein theory, such metrics will in general depend upon the choice of gauge. 

One way to resolve this ambiguity is to project the initial metric on SU (3) orthogonally to the 
orbits of U(1) 2 L x U(\) 2 R . The resulting Kaluza-Klein metric shall be discussed in Sec. 12.21 



2.1 The flag manifold SU(3)/U(1) 2 

A choice of metric on the coset SU(3)/U(1) 2 can give rise to different metrics on the bi-quoticnt, 
depending on the choice of section we make. In what follows, we shall illustrate this by choosing a 
natural metric on the flag manifold SU(3)/U(1) 2 that is Einstein-Kahler. 

There is a general construction showing that every quotient of a compact Lie group G by its 
maximal torus may be regarded as an Einstein-Kahler manifold. A physical application of this result 
would be to the modulus space of vacua of a Yang-Mills theory with Higgs in the adjoint. We shall 
describe the special case of G — SU(3), following the construction described in [17]. This makes use 
of the fact that 



S*[/(3)/E/(l) 2 = SX(3, 



B. 



(32) 



where B is the Borel subgroup of SL(3, 
the Iwasawa form 



In other words, we can express an SU(3) matrix U in 





( 


1 









(« 


u = 




-Z3 


1 












V 


~z 2 


Zl 


l ) 







(1 





yi 




1 


2/3 





1/ 



(33) 
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Substituting the expression for U given in (|16p . we find in particular that 
zi = —e~ 2iq tana; , 

Z2 = e~ 3lp ~ lq (e™ cos a; tany — sin a; secy tanz) , (34) 
Z3 = e~ 3lp+19 (e lt0 sin a; tany + cos a: secy tanz) . 

These expressions can be inverted to give the real coordinates in terms of the z a : 

\z 2 - zi z 3 | 2 



tan x — \zi\ . 



tan 2 y 



1 + M 2 ' 

2 



tan 2 z - lZ3 + ZlZ2 



1 + N 2 + \Z2-Z1 z 3 \ 2 
(z 3 tan x — z\ Z2) tan z 
(23 + z\ z-i tan a:) sin y ' 



with p and q then obtained using 



2 1 _ g-4ig £2_ _ e -6ip-2ig+2ito 

Zl ' Z2 



(35) 



(36) 



As discussed in [17], the z a can be viewed as complex holomorphic coordinates on the flag 
manifold. The Kahler function is given by 

K = log(l + |z 2 | 2 + |z 3 | 2 ) + log(l + |zi| 2 + \Z2 - z x z 3 | 2 ) • (37) 

It is easy to check that the Kahler metric, given by 



d 2 K 

ds 2 =g a pdz a dz\ g a p= gzag _- , (38) 



has a determinant given by 



X0 

which can therefore be written as 



det( ffQ a) = 2(1 + |z 2 | 2 + |z 3 | 2 )" 2 (1 + N 2 + \Z2 - zi z 3 \ 2 r 2 , (39) 



det(g a0 ) = 2e~ 2K . (40) 



Thus K satisfies the Monge- Ampere equation, implying that the Kahler metric g a p is Einstein. 
(Since R aB = d a dp\og(^/g).) 

Substituting (f34|) into (|38|) . one obtains the Einstein-Kahler metric on the flag manifold written 
in terms of the real coordinates (p, q, x, y, z, w). It is straightforward to verify directly that it satisfies 

Rij = igij . (41) 

In terms of the real coordinates, the Kahler function |3~7|) is given by 



lx — cos 2 x cos 2 z cos 4 y . (42) 
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The Einstein-Kahler metric (|38[) is invariant under the left action of SU(3>), and in particular, 
under the T 2 action generated by d/dp and d/dq. From this action corresponds to phasing the 
complex coordinates z a in such a way as to leave the Kahler function (|57|) invariant. It is possible, 
therefore, to perform a Kaluza-Klein reduction on the two angles p and q, to obtain a metric on the 
double coset U(1) 2 \SU(3)/U(1) 2 . The resulting metric is extremely complicated, and we shall not 
give it explicitly. However, the metric we obtain is not the same as the one discussed in Sec. 12.21 

This difference is connected with the fact that the Einstein-Kahler metric given by (|3"Tf and (|3"5|) 
is not the "round" metric on SU(3)/U{\) 2 , but rather, it is a particular member of a one-parameter 
family of homogeneous squashed metrics. (It corresponds to the only other member of the family, 
other than the round metric, that is Einstein.) The Einstein-Kahler metric constructed in (j3"5)) is 
given, in terms of the left-invariant one-forms a a defined in (|25p . by 

ds 2 = a\ + erf + cr 2 + cr 2 + 2 (cr 2 + er 2 ). (43) 

This illustrates the remarks we made previously about the ambiguity of measures on bi-quotients. 
The round metric on SU (3) corresponds to setting 

(a,/3, 1 ,S 1 ,S 2 ,d 3 ) = (1,1,1,1,1,0) (44) 

in (f29|) . Kaluza-Klein reduction with respect to d/dr and d/dt gives the round Einstein metric 
corresponding to a = f3 = 1 in (|30[) . The same metric on the flag-manifold quotient would also arise 
for general values of Si, S2, and $3, as long asa = /3 = 'y = l. If, on the other hand, 

(a, P, 7) = (1,1, 2), (45) 

for arbitrary Si, S2, and £3, we obtain the squashed Einstein-Kahler metric (|38j) on the flag manifold. 

This construction, while not providing us with a "simple" metric on the double quotient, has the 
virtue of being possible for any coset SU{N)/U{1) N ~ 1 ; we shall see in Sec. 16.21 that the case N = 6 
may be of relevance to neutrinos. 

2.2 Kaluza-Klein reduction of the bi-invariant metric 

Here, we start with the bi-invariant metric on SU(3), 

ds 2 = ±Tr dU dU^ = a 2 a . (46) 

In terms of the coordinates (p, q, r, t, x, y, z, w), it is given by 

ds 2 = 3dp 2 + dq 2 + 3dr 2 + dt 2 + f (3 cos 2y - l)dpdr + 3 cos 2 y (cos 2z dpdt + cos 2x dqdr) 
+^{cos2a; cos 2,2 (cos2y — 3) + 4sin2x sin 2,2 siny cosw} dqdt 

— sin 2 y (3dp — 3dr + cos 2x dq — cos 2z di)dw + 2 sin y sin w (sin 2z dtdx — sin 2x dqdz) 
+dx 2 + dy 2 + dz 2 + sin 2 y dw 2 + 2 sin y cos w dxdz . (47) 
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As expected, this metric on SU (3) does not depend on p, q, r and t, which are the arbitrary quark 
phases appearing in the Kobayashi-Maskawa matrix. 

It is perhaps worth remarking here that the first expression in (|46|) is well defined for any complex 
matrices U, unitary or not. For general complex matrices, it defines a flat metric on the space of 
matrix elements, which may be identified with C 9 = E 18 , the 18-dimensional Euclidean spaced 
SU(3) may be regarded as a real eight-dimensional submanifold of E 18 , defined by the nine real 
unitary constraints UW — 1 together with the one real unimodularity constraint det U = 1. The 
bi-invariant metric on SU (3) is the induced metric on this submanifold. 

One approach to placing a measure on mixing angles would be to give a uniform measure on the 
unconstrained mixing angles, and then to obtain a measure on the mixing angles by implementing 
the unitarity and unimodularity conditions. The left-invariant measure on SU(3) is unique up to a 
scale. Thus, any construction which respects SU (3) invariance will result in a measure which is a 
constant multiple of the Riemannian measure constructed from the bi-invariant metric. 

Writing (|47p in the standard Kaluza-Klein form, 

ds 2 = hij(x) (dy l + A i ^ l (x)dxi")(dy j + Ai{x)dx u ) + g^(x)dx^ dx v , (48) 

where y % = (p, q, r, t) and = (x, y, z, w). The metric on the bi-quotient is then given by 

dS 2 = g f , v (x)dx f " dx u . (49) 

The metric (gSJ) is once again rather complicated, and we shall not present it explicitly since we 
really only wish to calculate the Riemannian measure 

A 4 = dxdydzdw . (50) 

Noting that det g = det h det g, and that 

det g — sin 2 2x sin 2 y sin 2 2z cos 6 y , (51) 
we find that, after extracting an unimportant overall constant factor, 

det g = sin 2 2x sin 2 2z sin 2 y cos 4 y/F , (52) 

where 

F = (sin 2 2x + sin 2 2z) sin 2 y + |(5 cos 2y - 3) sin 2 2x sin 2 2z 

+| sin 4a; sin 4z sin 3 y cos w + i (3 cos 2y — 5) sin 2 2x sin 2 2z sin 2 y cos 2 w . (53) 

Note that one can alternatively obtain the four-dimensional metric on the bi-quotient by means 
of a T 2 Kaluza-Klein reduction of the round flag-manifold metric 

ds 2 = of + a\ + o% + of + oi + a 2 - , (54) 

2 Obviously, for k X k matrices, C 9 is replaced by C k . 



13 



which differs from the one used in Sec. 12.11 

We have calculated the Riemannian metrics of the two different metrics on the bi-quotients, and 
confirmed that these two four-dimensional measures are indeed different. Later, we shall demonstrate 
the dependence of the mean-square value of the Jarlskog invariant on the choice of squashing. 

2.3 Squashed Kaluza-Klein metrics 

As we noted earlier, not only is the four-dimensional double-coset metric obtained by Kaluza-Klein 
reduction nonunique, but also the associated measure is nonunique. In Sec. 12.21 we constructed the 
measure that follows from the T Kaluza-Klein reduction of the bi-invariant SU(3) metric to four 
dimensions, or, equivalently, the T 2 Kaluza-Klein reduction of the round six-dimensional flag metric 
([54|) . Here, we present the more general result for the measure on the double coset that is obtained 
by Kaluza-Klein reducing a one-parameter family of squashed flag metrics. Specifically, we take as 
our starting point the flag metrics 

ds 2 = a\ + a\ + a\ + a 2 + [a\ + cr 2 ) . (55) 

After Kaluza-Klein reduction, we find that the determinant of the four-dimensional metric g^ v 
is given, after again extracting an unimportant overall constant factor, by 

det g — sin 2 2x sin 2 2z sin 2 y cos 4 y/F , (56) 

where 

F = (sin 2 2x + sin 2 2z) sin 2 y + i(5 cos2y - 3) sin 2 2x sin 2 2z + X(J3 - l) 2 sin 2 2x sin 2 2y sin 2 2z 



+09 -1) 



4 sin y sin z cos z 



sin 2 2x | cos 2 z cos 2 2z — | cos 2z sin 2 2y — cos 2 y cos 4 z [cos 2 y — (3 + sin 2 y) sin 2 z] j 



sin Ax sin y sin 2z cos w 



cos 2z sin 2 y + (/3 — 1) (sin 2 y cos 2 z — sin 2 z) cos 2 z (57) 



+| sin 2 y sin 2 2x sin 2 2z cos 2 w 3 cos 2y - 5 - 4(/3 - l)(cos2z + sin 2 y) - 2((3 - l) 2 cos 2 y 

Note that this expression reduces to (f5"3")) if = 1, which is the special case of the reduction of the 
round flag metric. The nontrivial dependence of (|57p on the squashing parameter (3 shows that the 
measure yfg dxdydzdw on the double coset is also nontrivially dependent on the choice of squashing. 

2.4 Ozsvath-Schiicking metric 

The previous calculations give rise to rather complicated formulae. It is striking, therefore, that the 
metric obtained by Ozsvath and Schucking [9] is so much simpler. Their choice of section consists 
of simply setting p = q = r = t = in the metric (|47p . This results in the metric 

ds 2 = dx 2 + dy 2 + dz 2 + 2 sin y cos w dxdz + sin 2 ydw 2 . (58) 
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This metric is manifestly invariant under translating the coordinates x and z. Remarkably, there is 
one further commuting Killing vector. If one defines new coordinates u and v by 

(sin y cos w, sin y sin w, cos y) = (cos u, sin u cos v, sin u sin v) , (59) 

then the metric ([SHI takes the form 

ds 2 = du 2 + dx 2 + dz 2 + 2 cos u dxdz + sin 2 u dv 2 , (60) 

which has the three commuting Killing vectors d/dx, d/dz and d/dv. 

Geometrically, we can understand this if we note that the metric (|58p may recast as 

ds 2 = (dx + sin y cos w dz) 2 + (1 — sin 2 y cos 2 w) dz 2 + dy 2 + sin 2 y dw 2 , (61) 

which exhibits it as a T 2 fibration (having coordinates x and z) over a round hemisphere (having 
coordinates y and tu). (The colatitude y lies in the interval < y < \k ) The hemisphere can be 
embedded isometrically into Euclidean three-space as (siny cos w, sin y sin w, cosy). Equation (|59p 
then gives a different embedding such that while d/dw generates a rotation around the first axis, 
d/dv generates a rotation around the third axis. The projection along the third axis is given by 
cosy, while the projection along the first axis is given by cosu. 

Note that the extra Killing vector d/dv is purely local, since rotations about the first axis do 
not preserve the hemisphere. 

The metric (|S0"1) can be recast in the form 

ds 2 = sin 2 u dv 2 + (dz + cos u dx) 2 + du 2 + sin 2 u dv 2 , (62) 

which is locally of the form of a U(l) fibration (with coordinate v) over S 3 . 
The un-normalized measure is given by the remarkably simple formula 

fi = sin y(l — sin 2 y cos 2 w) 1 ^ 2 dxdydzdw , 

= sin 2 u dudvdxdz . (63) 

Despite its appealing simplicity, the Ozsvath-Schucking construction lacks a geometrical justi- 
fication and introduces a spurious U(l) 3 symmetry into the problem. A simpler example, which 
makes this clear, is provided by considering the lower-dimensional example of quotients of SU(2). 
The bi-invariant metric on ST/ (2) is 

ds 2 = (dip + cos 6 dcf)) 2 + d9 2 + sin 2 9 d<t> 2 , (64) 

where d/d<fi generates U(1)l and 8/ dip generates U(1)r. Projecting the metric orthogonally to the 
orbits of right translations, a la Kaluza-Klein, gives the round metric 

ds 2 = d9 2 + sin 2 6 d<p 2 (65) 

on S 2 . By contrast, simply setting dip = (the analog of the construction of Ozsvath and Schucking) 
instead gives the flat metric 

ds 2 = d9 2 + d<p 2 . (66) 
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The round metric (|65[) is invariant under £0(3). The flat metric (|66p appears to be invariant under 
the Euclidean group, with d/d9 and d/d<j> having the appearance of translations, but these are only 
local symmetries since is a periodic coordinate and 8 lies in an interval. 

The example of SU (2) also illustrates the difference between taking the flag-manifold measure 
and the Kaluza-Klein measure on a biquotient. The biquotient U(1)\SU(2)/U(1) = U(1)\S 2 is just 
an interval. Its metric becomes, after performing another Kaluza-Klein reduction of (|65p . 

ds 2 = d6 2 . (67) 

The measure would be dO, and not d9 sin 9 as obtained by integrating a function f{9) over the 
coordinate tf>. It is apparent from this simple example that there are inequivalent ways of calculating 
integrals of a function on a right quotient that is invariant under the left group action; namely, one 
can either reduce the metric to obtain a measure on the double quotient or take the measure on the 
single quotient and integrate out the left phases. 

3 Statistics of the Jarlskog Invariant J 

Expressed in terms of the coordinates (x, y, z, w), the Jarlskog invariant (|13p is given by 

J —\ sin 2x sin 2z sin y cos 2 y sin w . (68) 

The average of a function / on a space with metric is defined by 

^ = J fVgdxdydz dw 
J y/g dx dy dz dw 

The experimental value of the Jarlskog invariant J of the Kobayashi-Maskawa matrix is 

J = 3.08±g;Jf x 10" 5 , (70) 

which is very small compared with its maximum value 

J max = ^= w 0.0962 . (71) 

In the following subsect ions 1 3 . 1 1 1 o 1 3 . 5 [ we calculate the moments of J for the Kobayashi-Maskawa 
matrix, using the various measures we have introduced, and compare them with the experimental 
value. We will see that the average values one obtains are rather insensitive to the choice of measure. 

In subsection 13. U we start with 5[/(3)-invariant measures on the flag manifold SU(3)/U(1) 2 . In 
this case, as we have already noted, there is an unambiguous S'J7(3)-invariant measure. 

3.1 The flag-manifold measure 

In Sec. 11.21 we saw that the Kobayashi-Maskawa matrix is an element of the four-dimensional bi- 
quotient /7(l) 2 \S't7(3)/C/(l) 2 , which is however composed of two elements of U(1) 2 \SU(3). This 
left quotient is, just like the right quotient SU(3)/U(1) 2 , the flag manifold. 
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Since the Jarlskog invariant (fit)]) , or (f6"5| . is independent of all phasing angles, the averaging of 
its moments over the flag manifold will give the same results regardless of whether one constructs 
the manifold as the left quotient or the right quotient of SU(3) by U (l) 2 . This is convenient because 
we have already presented detailed results for the metrics on the right cosets SU(3)/U(1) 2 . 

A straightforward calculation shows that for the general class of 5J7(3)-invariant flag metrics 

yfg = 3a/?7 sin 2x sin 2z sin y cos 3 y . (72) 

Since an overall constant factor in the measure cancels out in the normalized averaging process, we 
see therefore that in contradistinction to the situation for the double coset U(1) 2 \SU(3)/U(1) 2 , the 
natural measure on the flag manifold is unique. 

One might think of using the Cartesian product of two flag manifolds for calculating the moments 
of the Jarlskog invariant for the Kobayashi-Maskawa matrix. In fact, one could equally well use 
SU(3) instead of the flag manifold since neither the bi-invariant measure nor J depend on the U(l) 2 
angles. The natural measure on SU(3) x SU(3), induced from is 

[i = N cti A <7 2 A (73 A cr 4 A ct 5 A cr 6 A 07 A cr 8 A a[ A a' 2 A <j' 3 A a\ A a' 5 A a' 6 A a' 7 A a' 8 . (73) 

Since it is only V = UU'^ that enters into the CP violating parameters, one could consider U 
and V as independent variables, i.e. write U' = V^U for some matrix V. Then the Maurer-Cartan 
form on the second SU (3) is 

\\ a a' a = U A dXJ' = UUU -rf(dVV*)U , (74) 

which gives a' a = a a — h a bTb, where 77, are right-invariant forms on SU(3) in terms of V coordinates 
and h a b only depends on the U coordinates. The measure (1731) . expressed in terms of V and U 
coordinates, is thus a product of a function of the U coordinates and the natural measure in V 
coordinates (left- and right-invariant forms on SU(3) give the same measure). Integration over the 
U coordinates then just gives an irrelevant constant, and one is left with the measure (|72[) on the 
space of V matrices. This justifies the use of ((T2"|) instead of the more complicated constructions 
obtained by reducing to the double quotient, and we will regard ()72|) as the most natural choice of 
measure on the parameter space. 

For the measure (|72")) the evaluation of the necessary integrals is very simple and we find that all 
odd powers of J average to zero, and 

(J 2 ) = — « 1.389 x 1CT 3 , (J 4 ) = — - — w 4.960 x 10" 5 . (75) 
x 1 720 . \ / 201600 v ; 

Thus we find that A J for the Jarlskog invariant is given by 

A J = —= « 0.0373 , (76) 
which is about three orders of magnitude larger than the experimental value (|70p . 
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3.2 The Kaluza-Klein measure from the bi-invariant metric 

For the metric on the biquotient discussed in Sec. 12.21 the expression for the measure is too com- 
plicated to allow us to perform the integrations analytically. Using numerical integration, we find 
that 

(J 2 ) w 1.1161 x 1CT 3 , (J 4 ) w 3.750 x 10~ 6 , (77) 
with the odd powers of J again averaging to zero. Thus we find 

A J rj y/Jj 2 ) ~ 3.341 x 10~ 2 , (78) 

which is very close to the previous result. 

Naively, one might have thought that since J is independent of all the U(l) phases, the results 
would be the same whether one averaged over the space U(1) 2 \SU (3)/t/(l) 2 , or else the flag manifold 
SU(3)/U(l) 2 . Of course we know that this is not in fact correct, since, as we have seen, the measure 
for the biquotient depends nontrivially on the squashing parameters a, /3, and 7 in (|30p while the 
measure on the single quotient does not. Nevertheless, it is interesting to compare the expressions 
(|75p for (J 2 ) and (J 4 ) with the ones we obtained in (177j) for the biquotient averaging. They are in 
fact quite similar, although the values are larger in (f75]j than in (I77[) . (We will see in the following 
subsection that, among the more general class of squashed biquotient measures, the values of (J 2 ) 
and (J 4 ) seem to be maximized by the round case (|77p.) 

3.3 The Kaluza-Klein measure from squashed metrics 

We can repeat the calculations of Sec. 13.21 using the measure given by (p^B"]) and ([5T|l for the one- 
parameter family of squashed Kaluza-Klein metrics. In view of the complexity of the measure, we 
must again resort to numerical integration. 

The case where the squashing parameter is chosen to be (3 — 2 is of particular interest, since 
this corresponds to the second Einstein metric on the flag manifold, i.e. the one associated with the 
Einstein-Kahler metric we discussed in Sec. 12.11 For this choice, we find 

(J 2 ) w 1.1012 x 1CT 3 , (J 4 ) w 3.678 x 10~ 6 , (79) 

with the odd powers of J averaging to zero. Thus we find 

A J w y/iJ 2 ) ~ 3.318 x 10 -2 , (80) 

This is smaller than the value of A J we obtained in (|78|) for the averaging over the Kaluza-Klein 
reduction of the bi-invariant metric, but only by about 0.7%. This does not bring it significantly 
closer to the experimental value of J, given in (|70[) . 

One might wonder whether, for some sufficiently large or small choice for the squashing parameter 
/?, it might be possible to obtain a result for AJ that was comparable with the experimentally 
observed value. In fact, it appears that (J 2 ) is a rather slowly varying function of /3. The value of 
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(J 2 ) appears to be maximized by the choice (3=1, and to fall off monotonically in both directions 
as /3 is taken to zero or to infinity. 

For example, if we choose (3 — \ we find 

(J 2 ) r* 1.103 x 1CT 3 , AJ« 3.321 x 1(T 2 , (81) 

while if we take (3 — > we find 

(J 2 ) ~ 8.097 x 1CP 4 , A J rj 2.846 x 10~ 2 . (82) 

Taking f3 = 1000, we find 

(J 2 ) w 4.298 x 10~ 4 , A J w 2.073 x 10~ 2 , (83) 

while for /3 = 10 6 we find 

(J 2 ) w 1.958 x 10~ 4 , AJ w 1.399 x 10~ 2 . (84) 
Even quite extreme values for the squashing parameter only bring about small reductions in A J. 

3.4 The Ozsvath-Schiicking measure 

Using the Ozsvath-Schiicking measure (|63|). we find that (J) = and 

(J 2 ) = 35 x 2~ 16 Ri 5.341 x 10~ 4 , (J 4 ) = 27027 x 2~ 34 w 1.573 x 10~ 6 , (85) 
and that the standard deviation is 



AJ 2 = V(</ 4 ) - {J 2 ) 2 = V22127 x 2~ 17 w 1.135 x 10~ 3 , (86) 

and 
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AJ- 7(J^ = A_« 2.311 x 10" 2 . (87) 
Again, the results are rather similar to the previous cases. 

3.5 The uniform measure 

Assuming a uniform distribution over the angles, and hence treating the double coset as a flat 
four-dimensional manifold so that the measure is simply fi = 1, would give 

(J)=0, (J 2 ) = t^tt ~ 4.883 x 10~ 4 , (J 4 ) = 189 x 2~ 27 w 1.408 x 10~ 6 . (88) 
2048 

Hence, this simplest possible choice gives 

AJ« 2.210 x 10~ 2 . (89) 
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4 Fine-tuning of J 



In the previous section we saw that different measures on the space of mixing angles all seem to 
lead to expectation values for J which are about three orders of magnitude larger than the observed 
value. The value for J that we observe hence appears to be finely tuned. In this section we shall 
do a closer, mainly numerical, analysis of the fine-tuning involved. We compare results obtained 
by taking the S'C7(3)-invariant and Kaluza-Klein measures, which seem natural from a geometric 
perspective, with a uniform distribution which is just the simplest possible choice. 

4.1 Probability distribution of J 

The observed value for the Jarlskog invariant J is 

J « 1(T 4 - 51 « e^ 10 - 39 . (90) 

In order to obtain a probability distribution for J we have used Mathematica to numerically compute 
integrals of the form 

J ^/g~dxdydzdw8(a-\J\)6(\J\-b)=P(b<\J\<a)- J ^gdxdydzdw (91) 

using Monte Carlo methods. The S'C/(3)-invariant flag-manifold measure and the Kaluza-Klein 
measure disfavor small values of J more strongly than a uniform distribution would. For example, 
we obtain 

fflagd J\ < 10~ 4 ) w 0.25% , P K k(| J\ < 1(T 4 ) w 0.44% . (92) 
Taking a uniform distribution y/g = 1, we get 

Punif(|J| < 10- 4 )«7%. (93) 



p(log(|/|» 




Figure 1: Probability distribution for log J using the SL r (3)-invariant flag measure, with fit top(| J|) oc | J\ x . 
The degree of fine-tuning required to reproduce a very small J is considerably higher if one uses the 
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measure induced by a £f (3)-invariant flag metric or the Kaluza-Klein metric, maybe contrary to 
what one might expect. Values of J close to its maximal value of w 0.0962 arc disfavored in 
both cases. Therefore we have used a logarithmic scale for \J\. 



p(log(|/|» 




Figure 2: Probability distribution for log |J| using the Kaluza-Klein measure, with fit to p(\ J|) oc | J\ x . 

In all three cases the numerical results for small |J| are well approximated by a power law of the form 
p(\ J|) = a ■ \ J\ X for the probability density of |J|. The logarithmic graphs show p(\og \ J\) oc | J| A+1 . 
For the 5'C/(3)-invariant flag measure (Fig. 1), the best fit to the data in the region below \ J\ = 10~ 23 
or log | J\ = —5.3 is 

A fl ag = -0.042(±0.006) , a flag = 18.1(±0.7) ; (94) 

for the Kaluza-Klein measure (Fig. 2) we fitted the data in the region below |J| = 10~ 2 - 7 or 
log|J| = —6.2 and obtained 

Akk = -0.097(±0.008) , a KK = 18.9(±1.0) ; (95) 

finally for the uniform measure (Fig. 3), the best fit to the data in the region below |J| = 10~ 3,4 

P(log(|/|» 




Figure 3: Probability distribution for log \ J\ using a uniform distribution, with fit to p(| J|) oc J\ x . 
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or log|J| = -7.8 is 

A unif = -0.500(±0.005) , a unif = 3.51(±0.15) . 



(96) 



4.2 Wolfenstein parametrization 

A different parametrization of the Kobayashi-Maskawa matrix which is frequently used was intro- 
duced by Wolfenstein and is based on the experimentally observed hierarchy 

y < x < z < 1 (97) 

in the mixing angles. One rewrites [18] 

sin z = A , sin x = AX 2 , sin ye~™ = AX 3 (p - vq) (98) 

and treats A as a small parameter while A, p, and r\ are supposed to be parameters of order unity. In 
the modern literature one also frequently uses p, fj instead of p and rj because then the combination p+ 
if) is independent of the phase convention in the Kobayashi-Maskawa matrix [19]. These parameters 
are defined by 



1-A 2 A 4 ^AV+jj l-A^ fj 

1 - A 2 (1 - A 2 X A p) 2 + A 4 AV ' ' V 1 - A 2 (1 - A 2 A 4 p) 2 + A 4 A 8 ?f ' [ ' 

The experimental values for A, A, p, fj are [19^ 

A = 0.2272 ± 0.0010, A = 0.818±g;g?£, p = 0.221±g;g||, fj = 0.340^;^ . (100) 

One viewpoint on the Wolfenstein parametrization is that it is adapted to the values for the 
Kobayashi-Maskawa matrix entries that we observe and has no deeper significance; but often the 
viewpoint is expressed that this parametrization expresses some kind of "natural hierarchy" in the 
mixing angles coming from physics beyond the standard model (see e.g. [20]). Treating the other 
parameters as "naturally of order unity" reduces our calculations to a one-dimensional problem as 
everything is only expanded in terms of A. We find that the 5?7(3)-invariant measure on the flag 
manifold is now, to leading order in A, 

d(x,y,z,w) 



d(X,A,p,ri) 
and the Jarlskog invariant J is 



y/g oc A 3 X n (l + A 2 + 0(A 4 )) , (101) 



j = A 2 yX%l - yl 2 A 4 ) (l - A 2 - 2v4 2 pA 4 - A\^ + (p- 2)p)X e + A\n 2 + P 2 )A 8 ) = j4 2- A 6 +0(A io ) 
(1- A 2 )(l-2A 2 pA 4 + A 4 (?7 2 +p 2 )A 8 ) 2 V 

(102) 

Inverting this expression to leading order gives the probability distribution for J 

p{J)(X ^f( 1+ {^) 1/3 + 0{j2/3) ) ' (1 ° 3) 



3 Note that only even powers of A appear in all expansions, so that it is A 0.05 which is the small parameter. 
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which is incompatible with the numerical results. Trying to improve this approximate result by 
letting A, p and fj take all possible values leads to inconsistencies since the expansion in powers of J 
contains poles of arbitrary order in A. From our present viewpoint, where no mechanism for fixing 
these parameters close to one is known, the Wolfenstein parametrization seems rather misleading 
when discussing geometric probability. 



5 Quark Mass Matrices and Gaussian Weighting Functions 

In the previous sections we have focussed on U(l) 2 \SU{3)/U(l) 2 , the space of Kobayashi-Maskawa 
matrices, as the space of CP violating parameters. Since SU{3) is compact, this space has finite 
volume for a natural measure. But the Kobayashi-Maskawa matrix is derived from the Hermitian 
quark mass matrices, which could be viewed as more fundamental and more directly determined 
by physics beyond the standard model. In this section, we try to obtain statistics of the Jarlskog 
invariant J from a random distribution on the space of 3 x 3 Hermitian matrices. 

5.1 Distributions on Hermitian matrices 

We follow Sec. 11.21 and write the quark mass matrices as 

UmU^ = diag(m u , m c , m t ) , U'm'U = diag(mrf, m Sl m^) . (104) 

where U and U' should be thought of as elements of U{1) 2 \SU{3). Following [12] we normalize the 
mass matrices by dividing by mass scales A and A' (often taken to be the top and bottom quark 
mass, respectively) which may be chosen for convenience: 

M = U^DU , M' = U A D'U' . (105) 

The matrices D and D' are now dimensionless quantities, and it is clear that U and U' are only 
defined up to left multiplication by elements of U(l) 2 . We consider A and A' as arbitrary mass 
scales, and so we will allow arbitrary eigenvalues for both matrices, instead of fixing one of them to 
be one. 

A natural measure on the space of Hermitian matrices is induced by the metric 

ds 2 = Tr(dM ■ dM) = Tr {dD ■ dD) + 2Tr {{dU t/ f D) 2 - (dU U^f D 2 } (106) 
which is invariant under conjugation under U(3). If we define right-invariant one- forms r a by 

dUU^ =iA a r Q , (107) 



this becomeqj [with D = diag(Z?i, D 2 , D 3 )] 

ds 2 = Tr {dD ■ dD) - 2r a r b Tr (A a [D, X b ]D) 

= dD 2 + dD 2 + dD 2 + 2 {(£>! - D 2 ) 2 {t 2 + t 2 ) + {D l - D z ) 2 {tI + t 2 ) 

+ {D 2 -D 3 ) 2 {r 2 +r 2 )} . (108) 

4 Compare with the corresponding result for real matrices given in [21] 
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The corresponding volume form is 



(Di - D 2 ) 2 {D 1 - D 3 f{D 2 - D 3 f dDi A dD 2 A dD 3 A n A r 2 A r 4 A r 5 A r 6 A tj . (109) 

As explained above, the measure on the coset U(1) 2 \SU(3) is unique and equal to the measure 
induced from the bi- invariant metric on SU(3). We obtain a Riemannian measure 

DM := (L>! - D 2 ) 2 (£>i - D 3 f{D 2 - L> 3 ) 2 sin2x cos 3 y siny sin2zdD 1 dD 2 dD 3 dx dy dz dw dr dt 

(110) 

on the space of Hermitian 3x3 matrices. The coordinates (x,y, z,w,r,t) on U(1) 2 \SU(3) were 
introduced in Sec. and we allow arbitrary eigenvalues. (In a fermionic mass term, the sign of the 
mass has no physical significance, since it can be reversed by multiplying the spinor fields by 7 5 ; 
only m 2 enters in physical quantities.) 

From the expressions (|105[) , it is apparent that each Hermitian matrix with three distinct eigen- 
values is associated with six different elements of R 3 x U(1) 2 \SU(3), related by the action of the 
discrete group 63: 

M = U^DU = (U^ P~ 1 )PDP~ 1 (PU) =: U^DU , P G 63 , (111) 

where S3 is the symmetric group of degree 3 (the dihedral group of order 6, sometimes denoted by 
D 3 or Dq) which permutes the canonical basis vectors of M 3 . The set of matrices with coinciding 
eigenvalues has zero measure and hence can be ignored in the present discussion. 

Thus we need to consider the space E 3 x {U{1) 2 x & 3 )\SU(3) instead, restricting the coordinates 
on the flag manifold to an appropriate range to pick one of the six matrices related by the 63 action. 
We can use the fact that the 63 action permutes the rows of an SU(3) matrix to demand that the 
elements of the third column (see Sec. [2|) satisfy the relation 

I siny| < I sin 2; cos ?/| < | cos xcosy\ , (112) 

which restricts the coordinates x and y to 

7T 

< y < arctan(sina;) , < x < — . (113) 

Using the natural measure on the flag manifold, we see that this region has precisely one-sixth 
of the total volume of the flag manifold: 

7r/2 7t/4 arctan(sinx) 

J dz J dx J dy sin 2x cos 3 y sin y sin 2z 

— °- °- = -. (114) 

7r/2 tt/2 tt/2 6 

J dz J dx J dy sin 2x cos 3 y sin y sin 2z 
000 

An integral over R 3 with the given measure diverges. We could introduce a cutoff for the quark 
masses, but then any expectation values for quark masses would strongly contradict observation, as 
there is no way to explain the observed mass hierarchy. 
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We therefore choose to introduce a weighting function in the measure which decays sufficiently 
fast for large positive or negative eigenvalues and is able to reproduce the known hierarchy. The 
simplest assumption is to take a weighting function of the form 

/ (Tt(M 2 A)) / (Tr((Af ) 2 A')) , (115) 

where A and A' are Hermitian and positive definite, and we shall further assume [A, A'] = 0. 
By a redefinition of M and M' by unitary conjugation by the same unitary matrix, which leaves 
J invariant, one can simultaneously diagonalizc A and A' . For simplicity and ease of technical 
calculations, we shall choose the function / in (|115[) to be a decaying exponential so M and M' are 
governed by Gaussian distributions. Our proposal is to fit the diagonal matrices A and A' to the 
observed quark masses and use the resulting probability distribution for statistics of J. 
An integral of a quantity such as J 2 becomejfl 



(J 2 ) - NjDMDM'e- T ^ MA ^^ M 'r A ^J 2 (M,M') (116) 
= N [ dDdD 1 [ DUDU'e- T < D2uAU '~>~ Tl( ^ 2u ' A ' u '^J 2 (U,U'). 

Here DU and DU' are the measures on (U(l) 2 x & S )\SU(3) and dD := (D x - D 2 ) 2 (D 1 ~ 
Dz) 2 {D 2 — D3) 2 dDi dD 2 dDs etc., and the normalization factor N is defined by 



1 

N 



dDdD' DUD tf e -^&VArf)--mn?)*u'A>u'i)_ (n7) 

J((c/(i) 2 xe 3 )\sc/(3)) 2 



From J = 3m(Vii V 22 V* 2 V 2 \) and V = UU'\ we have 

3 

J(u,u')= Y, M u ^U2 b u^u; d u , * la u'; b u' 2c u' ld ). (us) 

a,b,c,d— 1 

At this point, it is perhaps instructive to note that setting A equal to the identity would split 
the integral (|116p into a product of an integral over the eigenvalues which just gives a constant and 
an integral of J 2 over ((U(l) 2 x &3)\SU(3)) 2 . Since all even powers of J are invariant under the 
S3 action on U and V , this can be replaced by an integral over (U(1) 2 \SU(3)) 2 if averages are 
concerned. By the arguments presented in Sec. 13. 1[ a change of coordinates reduces this to a single 
integration over a flag manifold, and one recovers the results of Sec. 13.11 for expectation values of 
powers of J. 

The introduction of more general diagonal matrices A and A' means that the invariance of the 
measure DAI DM' under separate conjugation of M and M' by arbitrary elements of U (3), i.e. under 
the action of U(3) x i/(3), is broken down to the action of the diagonal subgroup U(l) 2 x U(l) 2 
which commutes with A and A'. We find that this symmetry breaking is necessary to obtain a 
distribution that reproduces different expectation values for squared quark masses. 

It should be clear from (|116p that multiplying A (or A') by a constant is the same as rescaling 
the eigenvalues Di (or D-) and so amounts to a rescaling of A (or A'). We can therefore, without 



3 All odd powers of J again have expectation value zero. 
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any loss of generality, choose 



.4 



(\ \ 

1/nl 
\0 l/fil) 



A' 



A o o \ 
o o 

V° V/*3/ 



(119) 



where fj, c , fi u , fj, s , and /x^ are dimensionless parameters that we are free to choose so as to reproduce 
the observed quark masses as expectation values. (In the case of an exponential exp(— Tr(D 2 A)), 
these would of course be equal to the respective quark masses, expressed in units where A = nit and 
A' = rrib.) Because of experimental uncertainties in the up and quark masses, one can modify this 
distribution to reproduce different values for these masses. 

It seems practically impossible to evaluate the integral (|116j) , as the expression for J in terms of 
coordinates on ((C/(l) 2 x &3)\SU(3)) 2 is too complicated to be given explicitly. However, since 

Tr(D 2 UAtf) = J2D 2 a J2 A o\Uac\ 2 =: £ ^ , Tr((D') 2 U' A' U'^) =: J^jVa (120) 

a c a a 

with 

£i = A\ cos 2 y cos 2 z+A 2 cos 2 y sin 2 z+A$ sin 2 y , = A\ cos 2 y' cos 2 z'+A' 2 cos 2 y' sin 2 z'+A 3 sin 2 y' , 

(121) 

and we assume A3 3> 1 and A 3 ^> 1, the integrand is negligibly small unless y w and y' w 0. We 
use this to approximate the integrals over j/ and y': 

arctan(sin x) arctan(sin x') 

J dy J dy' cos 3 ysiny cos 3 y' smy'e- T ^ D2uAu ^- T ^ 2u ' A ' u ' f Kj 2 (U,U') 


arctan(sin x) arctan(sin x') 

dy j dy'yy'e-^ 2 -^ 2 ( e -^UAV*)-Tr«D>)*u>A>u'*) ja^y^ ■ 



y=y'=o 
(122) 



It turns out that this is independent of u> and iu'. Constant prefactors such as l/AA^A'^ appearing 
in both numerator and denominator can be dropped, and so we have 



(J 2 ) 



where 



/ R , dDdD' Jd 4 xjd 4 x' Sin2xsin2zsin2x'sin2z' ( e -i*(&UArf)-Tt«D') a u'A'u'i) J*(U,U') 



\y=y'=0 



/ R6 dD dD 1 J d 4 x J d 4 x> sin 2x sin 2z sin 2x> sin 2z' ( e -Tr(^ t /AC/t)„ T r((D') 2 c/'A'c/'t )) | ^ 

(123) 



tt/4 tt/2 2tt 2tt 

<i 4 :r = I dx I dz I dr I dt 



(124) 



0000 



and similarly for J d 4 x' . 
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Now we can integrate over both copies of M 3 in (|123|) . using 

oo oo oo 

J dD 1 J dD 2 J dD 3 {D x - D 2 f{Di- D s f{D 2 - D^fe-^-^-^ 

— oo — oo — oo 

3tt 3 / 2 

= ,.5/2,5/2.5/2 fe'fe + 6) + c1(6 + 6) + £1(6 + 6) - 2666) • (125) 



The explicit expression for J at y = y' = is 
I _ 1 

ly=v'=0 4 



^OL^^o = \ s 2 X s 2x , {c 2 z ,s 3 z s z , sin(3f + t) + clc z , S \, sin(3r - t) 



-<? z s z s z \<? z , [sin(3f + i) + sin(3f - 3t)] - s^, sin(3f + <)) 
+c z c zl s 2 z {c 2 zl sin(3f - £) - s 2 , [ s j n (3f + 3<) + s in(3r - *)])} (126) 

where = sin a;, cy = cos z', etc., f — r ~ r', and i — t — t' . Integrating (|126p over r, r', t, and i' 
indeed gives zero, which is why we choose to use J 2 . 

5.2 Results and dependence on quark masses 

We need to determine the parameters appearing in the matrices A and A' in (|119[) . We first observe 
that expectation values for squared mass matrices take the relatively simple form 

7r/4 tt/2 



J K 3 dD D\ f dx J dz sin 2x sm2z (e^ 13 * UAU '- 

n n \ 



ly=Q 



(D'() » r r : : • (127) 

1 tt/4 it/2 

J R3 dD J dx J dz sm2xsin2z(e-^( D2UAU1 )) \ 



The denominator is explicitly 

7r/4 tt/2 



I D := / dz / dz sin2zsin2z (6 2 (6 + 6) + 6 2 (6 + 6) + 6'(6 + 6) - 2666) , 

(128) 



where 



6 = A\ cos 2 z + A 2 sin 2 z , 6 = ^i c °s 2 x sin 2 z + A 2 cos 2 a; cos 2 z + A 3 sin 2 x , 

6 = Ai sin 2 x sin 2 z + A 2 sin 2 a; cos 2 z + A 3 cos 2 x , (129) 

with A3 ^ A 2 ^>> Ai. We notice that all 6 are nonzero for all values of x and z. Furthermore, the 
integral is dominated by very small x and z (we cannot have x = and we can approximate Id 
well by only keeping the terms of leading order in x and z in the trigonometric functions, and 

6 2 (6 + 6) + 6 2 (6 + 6) + 6 2 (6 + 6) - 2666 « A\ x 2 + A\A 2 , (130) 
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which are the leading terms (as we shall see, the first of these is effectively also of order A\A 2 )\ 

7r/4 7t/2 

3tt 3 / 2 



Id « ^ J dx J dz4xz(Alx 2 +AjA 2 )(A 1+ A 2 z 2 )- 5 / 2 (A 2 + A 3 x 2 )- 5 / 2 A- 5/2 



OO CO 

q 3/2 /■ /■ 

« dX dZ (A 3 3 X + A§^ 2 )(^i + A 2 Z)- 5 / 2 (A 2 + yl 3 X)- 5 / 2 ^ 5/2 

3tt 3 / 2 



o o 

OO 



■ /rfX(A 3 X + ^ 2 )(A 2 + A 3 X)- 5 / 2 V /2 -— i- 
8 J 3 Ay A 2 

o 1 



^2 ^3 _ .3/2 .1/2 I _ _ .3/2 .3/2 .3/2 ' V 1 - 31 ; 



4Af%V 3 2 3 3AfA 2 /V 2Af 2 /lf4 /r 
Similarly, we find 

15^ 3 / 2 



OO OO 



WD « 1 -^JdxJdZ(AlX + A 2 3 A 2 )(A 1 +A 2 zr 7 / 2 (A 2 + A 3 X)- 5 / 2 A- 5/2 
o o 

3/2 °° 

2 /dX(A 3 X + A 2 yl 2 )(yl 2 + A 3 X)- 5 / 2 A- 5 / 2 .^ 7I 
J bA 2 A-y 
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o 

37T 3 / 2 



hence 



(132) 



<^>«2^- (133) 



Redoing the same calculation for D 2 and D 3 gives 

(^-sb- <D|> "5i- (134) 

There is a relative factor of 3 which has to be taken into account when determining A and A' . 

Because of the dependence of masses on the energy scale in quantum field theory, described by 
the renormalization group, there is some ambiguity in what is meant by the "quark masses" we want 
to reproduce. Following [22], for example, we take all the quark masses evolved to the scale of the 
Z boson mass. These are given in [23]: 

(m u , m c , m t ) = (1.27±g;|^ MeV, 0.619 ± 0.084 GeV, 171.7 ± 3.0 GeV) ; 

(m d , m s , m h ) = (2.90ti; 2 g McV i 55 ±i5 McV > 2 - 89 ± °- 09 GeV) . (135) 

We use the central values 

(m u ,m c ,m t ) := (1.27 MeV, 0.619 GeV, 171.7 GeV) ; (m d , m s , m b ) := (2.9 MeV, 55 MeV, 2.89 GeV) 

(136) 
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The mass scales A and A' are fixed by setting (D\) = (m t /A) 2 and ((D[) 2 ) = (r7i fc /A') 2 . By 
comparing the results obtained by numerical integration with the values we want to reproduce, we 
can then fix the parameters fi c , fi u , /i s and /id- 

In the case of the positively charged top, charm and up quarks, which exhibit a more extreme 
quark mass hierarchy, we find that numerical calculations (using Mathematica) reproduce the results 
we have obtained analytically very well (see Table [TJ. For the negatively charged quarks, we find 
numerically that we have to use relative factors different from 3 to reproduce the observed masses. 
Comparing the numerical results with (|I36|) . we fix the parameters appearing in A and A' to 

^ = 3 ps 3.90 x ICT 5 , M« = 3 (—} « 1.64 x ICT 10 , 

\m t J \ m t J 

A = \[— V-5.43XH)- 4 , Md = T~ (^-\ 2 Ri 2.42 x 10~ 6 . (137) 
2 \m b ) 5 \m b J 

As a brief side remark, we see that the dominant contributions to these integrals come from the 
regions 




(138) 



and these values are all small compared to one. We can therefore give rough estimates for magnitudes 
of individual elements of the Kobayashi-Maskawa matrix. 

In the standard convention the ordering of the quark families is (u, c, t) and not (t, c, u) as used 
in (1119p . which means that in our parametrization, 

1(^)131 = 1^1, \(UU' i ) 12 \=\V ts \, \(UU'%\=\V cd \. (139) 

Since all of the numbers in (|138p are small, we only keep leading terms in the angles on U and U': 

\(UU^) 13 \ w \x'{z' -z)- e iw 'y' + . . . | w x'z' « W «2x 1(T 3 , 

\(UU'%\ «z'« Ms «0.02, |(^ /t ) 33 |wa/«^»0.07. (140) 

Ms 

Experimental values are [19] 

\V td \ = (8A4t° il) x 10~ 3 , \Vu\ = (41.61±g : ^) x 10~ 3 , \V cd \ = 0.2271^;°^ • (141) 

Our rough estimates reproduce the right ordering of the three parameters and are accurate to 
factors of order a few. A more careful analysis would involve computing expectation values for these 
parameters in the distribution we have assumed. 

We return to the task of computing the expectation value of J 2 . In order to obtain an analytical 
expression, we use the fact that the main contribution to the integral (|123j) will come from small z 
to only take the term in (|126[) that is nonzero at z = 0. Averaging over r, t, r',and t' gives a factor 
of 1/2, as one might have expected, and therefore we use 

J s 2 mall z := — sin 2 x cos 2 x sin 2 x' cos 2 x' cos 2 z' sin 4 z' (142) 
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for our calculations. The numerator of (|123[) is the product (using again that only small z contributes) 

tt/2 tt/2 

9n 3 f , 2z f , , sin 2z' cos 2 z' sin 4 z' 



x 



dz — : Trm^ X / dz' 



32 J (A 1 +A 2 z 2 f/ 2 J (Ai cos 2 z' + A 2 sin 2 z') 5 / 2 



tt/4 





tt/4 



, sin 2x sin 2 x cos 2 x(A2 cos 2 x sin 2 x + A 2 -A 2 (cos 6 x + 2 cos 2 x sin 2 x + sin 6 x)) 

<ix — — — — (143) 

(A 2 cos 2 x + A 3 sin 2 x) 5 / 2 (A 3 cos 2 x + A 2 sin 2 x) 5 / 2 



, sin2x' sin 2 x' cos 2 x'^A^) 3 cos 2 x' sin 2 a;' + (A' 3 ) 2 A' 2 (cos & x' + 2 cos 2 x' sin 2 x' + sin 6 x')) 
(A' 2 cos 2 x' + A' 3 sin 2 a;') 5 / 2 (A' 3 cos 2 x' + 4 sin 2 x') 5/2 

The first two factors are 2/(3A^ 2 A 2 ) and 4/ (3y/~A^(y/~A^ + V^A^) 4 ), respectively; for the other two 
(which are identical) we change variables to X = cos 2 x to obtain 
l 



f X(l - X)(A 3 X(1 -X) + A 2 A 2 (X 2 -X + l)) 
J (A 2 X + 4,(1 - X)f/ 2 (A 3 X + A 2 (l - X)f/ 2 ~ A 2 I 

1/2 V 



arctan i# -2, 
2 V 4 



(144) 



where we are dropping corrections of order Putting everything together, we obtain 



(■Call 2 > « (£)^yg2 / arctan f^f _ /4A7\ / arctan l£L_ I ^2(1, , , 



A' 3 4A' ? 

V4^(V^+V^) 4 V V4 ^ 2 v A3 A V 4 ^ 

^5 m « m ^ m ' / m c 2m u \ ( rYm s [32 m d 

| arctan I arctan - ' 



2 , /2iwV V 2to « m c J \ \ 32m d V 5 m 

m c m 2 (l + v /|^j \ 



where the numerical factors appearing in the last line come from the different factors chosen in (|137p . 
Note that the top quark mass does not appear in this approximate result. 

For numerical calculations we use both the simplified expression J 2 mal i z and the expression for 
J given in (|126p . We find that for the first quantity, the numerically evaluated expectation value 
( J 2 mall z ) is about 7/6 of (|145|) , and the numerical result for (J 2 ) (taken at y = y' = 0) is 

(J 2 ) ps 5.28 x 1CT 9 , (146) 

which gives 



A J = y/JJ 2 ) w 7.27 x 1CT 5 (147) 

which is much closer to the observed value than any of the previously obtained results. Assuming 
a Gaussian distribution for J which is peaked at zero, the probability of finding a small J, in the 
sense of Sec. |H is now 

P mass (|J| < HT 4 ) -83%, (148) 
whereas the probability of finding a J which is even smaller than the observed value is 

Pmass(|J| <3x 10" 5 )«32%. (149) 

The observed value for J can no longer be viewed as being finely tuned if the distribution used in 
our calculations is assumed. 
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Table 1: Analytical and numerical results for integrals of interest. 



Quantity 


Id (over x, z) 


Id(DI) 


Id{DI) 


Jd(DI) 


Analytical result 


1.43 x 10 


r\ -i a 1 a — 51 

2.14 x 10 


2.78 x 10 


1.17 x 10 J 


Numerical result 


1.43 x 10 


2.14 x 10 


2.78 x 10 


1.18 x 10~ 31 


Quantity 


(over x', z') 


^((•D'l) 2 ) 


I' D {{ D 2?) 


-r 1 t / I ""« / \ O \ 

^((^3) > 


Analytical result 


1.32 x 10 J 


1.99 x 10 a 


3.60 x 10 ' 


1.60 x 10 a 


Numerical result 


i or* -i /a — 1 1 

1.32 x 10 J 


1.98 x 10 J 


7.28 x 10 ' 


1.99 x 10 J 


Quantity 


Id (over x, z, ar , 2 J 


/ t2 \ 
Wsmall 2/ 


/ t2\ 
W ) 




Analytical result 


1.89 x 10~ 34 


3.22 x 10~ 9 






Numerical result 


1.88 x 10~ 34 


3.73 x 10- 9 


5.28 x 10" 9 





To test the sensitivity of our results to changes in the parameters, we take values at the upper 
or lower limit in (|135jl and try to find the highest and lowest values for (J 2 ). We find that setting 

(m u ,m c , m t ) := (0.85 MeV, 0.535 GeV, 174.7 GeV) ; (m d , m s , m b ) := (1.71 MeV, 40 MeV, 2.98 GeV) 

(150) 

gives 

(J 2 ) w 1.86 x 10~ 9 (151) 

and 

A J = VT^>" ~ 4.31 x 10~ 5 , (152) 

whereas setting 

(m», m c , m t ) := (1.77 MeV, 0.535 GeV, 168.7 GeV) ; {m d , m s , m b ) := (4.14 MeV, 71 MeV, 2.8 GeV) 

(153) 

gives 

(J 2 ) w 1.52 x 10~ 8 (154) 

and 

AJ = v 7 ^ 2 !- 1-23 x 10~ 4 . (155) 

Since the former choice makes J appear more typical, these results are perhaps an indication that 
the correct values for the (up, down, and strange) quark masses are probably closer to the lower than 
to the upper bounds given in Q135p . Also, even the greatest possible value for A J is significantly 
lower than any of the values obtained in previous sections. 

In this section, we have established that assuming the observed hierarchy in quark masses in a 
Gaussian distribution over the space of mass matrices gives expectation values for J 2 which are small 
enough to regard the observed value as "natural" and not finely tuned. This statistical observation 
seems to open up the possibility that the same mechanism that is responsible for the apparently 
unlikely hierarchy in quark masses might also explain why the observed value for \J\ is so small. 
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A more detailed analysis including a probability density for \J\ for this distribution is left to 
future work, since the numerical methods used here do not give sufficiently accurate results. 

6 Extension to Neutrinos 

In this section we review the case of neutrino masses, highlighting the difference between Majorana 
and Dirac masses and commenting on some recently made suggestions in the literature that one 
could distinguish the two cases by gravitational effects. In contrast to the physical situation which 
is at present rather unclear, the mathematical problem of obtaining a measure on the space of mixing 
matrices is in this case simpler, since one considers the flag manifold U(l) 2 \SU(3). 

6.1 Neutrino masses 

The spectrum of neutrinos and their masses and their nature, Majorana or Dirac, is currently not 
well known. Since not all of the material may be familiar to all readers we shall review some basic 
facts about Dirac and Majorana masses in a framework which is sufficiently general to encompass 
all likely possibilities. We find it helpful to use a Majorana notation, but no loss of generality 
thereby results since, if one starts with complex Weyl notation, one may always take real and 
imaginary parts. Alternatively, given a treatment in terms of Majorana spinors, one may always 
transcribe it into Weyl notation. In order to simplify the analysis we depart from common practice in 
phenomenological particle physics and adopt the spacetime signature ( — h ++) for all spinors. This 
has the advantage that all gamma matrices may be taken to be real, * denotes transpose. C = — C* 
is the charge conjugation matrix and j 5 — — (7 s ) 4 = 7°7 1 7 2 7 3 so that (7 s ) 2 = —I. If required, a 
concrete representation is given by 
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(156) 



In this representation we may take C — 70, and it is often useful to note that 7 is antisymmetric 
while 7 X ,7 2 and j 3 are symmetric. 

The mass matrix of a set of fcrmions is defined by imagining setting to zero that part of the 
effective or large distance Lagrangian containing couplings to all gauge interactions except gravity. 
In fact this is how information about neutrino masses is obtained. One observes mixing as they pass 
from the sun to the earth and upper bounds on their masses have been obtained by observing their 
arrival times from distant supernova 1987a. 
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The most general effective Lorentz-invariant Lagrangian for a system of k free Majorana fermions 
if) 1 , i = 1, 2, . . . , k, is 

L = ^ t GpV-^ t C(m 1 +m 27 5 )V, (157) 

where should be thought of as a fc-dimensional column vector all of whose entries are four com- 
ponent real Majorana spinors, and mi and rri2 are real symmetric k x k matrices. Note that at this 
stage k may be even or odd. We have made use of the fact that we may diagonalize the kinetic term 
using GL(k,R) transformations. We are still allowed SO(k) C GL(k,R) transformations 

ip^Otp, 0*0 = 1, (158) 

where 

O = cxp (wij) , u>ij = -ujji . (159) 
The kinetic term, but not the mass term is also invariant under chiral rotations 

ip^Pip, (160) 

P = exp (fjj7 5 ) , Vij = Vji . (161) 

Combining these two sets of transformations we see that the kinetic term, but not the mass term is 
in fact invariant under the action of U(k), i.e. under 

ip^Sip, (162) 

S = cxp (w I3 + . (163) 

The U (k) invariance is perhaps more obvious if one uses a Weyl basis. Since 

( 7 5 ) 2 = -l, (164) 

one may regard 7 5 as providing a complex structure on the space of 4fc real dimensional Majorana 
spinors, converting it to the 2k complex dimensional space of positive chirality Weyl spinors for 
which 

7 5 =i. (165) 
Clearly S then becomes the exponential of the k x k anti-Hermitian matrix 

ujij + ii/ij . (166) 

Thus 

55 f = 1 . (167) 
The mass matrix is then a complex symmetric matrix 

in = vn\ 7 im 2 , (168) 

and under a U(k) transformation 



m 



S^mS. (169) 
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At this point we invoke the result of Zumino [24] that S may be chosen to render the matrix m 
diagonal with real non- negative entries m». 

In the general case, all the masses rm are distinct. They are then said to be of Majorana type. 
However, it may happen that two masses, mi and mi say, coincide. One may then combine ip 1 ^ 2 
into a Dirac spinor. One then has the case of a Dirac mass. For quarks and charged leptons all 
masses are of Dirac type. For neutrinos, however, it is not yet known of what type they are, nor 
indeed how many. A simple assumption is that k = 6, with three having very heavy masses and 
three having very light masses. This corresponds to the so-called "seesaw mechanism" . Of course 
the very light neutrinos may be combined into three Weyl neutrinos and are taken to be massless in 
the standard model. 

From the analysis above it follows that in the general case when all masses m, are distinct, there 
is a unique basis for the neutrino states, determined by their inertial motion. If, however, two or 
more masses coincide, then the basis becomes ambiguous up to rotations of the components with 
equal masses. In the case that k is even and there are [|] distinct pairs of coincident masses the 
basis is arbitrary up to the action of T^] = U(l)^. 

Any mixing matrix taking one to a basis which is preferred from the point of their nongravita- 
tional gauge interactions will be ambiguous to the extent that the inertial basis and the gauge basis 
are ambiguous. For that reason, in general, a mixing matrix belongs to a double quotient. 

6.1.1 The Universality of free fall 

In our discussion above we have referred to the interactions of neutrinos with gravitational fields. 
Of course neutrinos observed to be coming from the sun, or the supernova 1987a are traveling so 
fast that the effects of gravity on them are negligible. However, it has been suggested that it is in 
principle possible to distinguish Majorana from Dirac masses by their behavior in the gravitational 
fields of rotating objects [25-27]. Our analysis above shows that unless there are gauge interactions 
such as might correspond to neutrino magnetic or electric dipole moments this is not so, as long as 
the coupling to gravity is "minimal." If so one simply uses for V the standard Levi-Civita covariant 
derivative acting on spinors. 

Assuming that the mass matrix m is independent of position, we may take it to be everywhere 
real and diagonal. Thus each component -0* of the inertial basis propagates independently. One may 
iterate the Dirac equation and use the cyclic Bianchi identity in a curved space to get (reinstating 
powers of h) 

- ft 2 V V + ^h 2 R^ + m 4 V = . (170) 
If the effects of curvature are negligible on the scale of the Compton wavelength, 

<™> 

the second term may be dropped and one obtains the Klein-Gordon equation for each component. 
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As is well known, there is no "gyro-magnetic" coupling between the spin and the Ricci or Ricmann 
tensors [28]. To proceed, one may pass to a Liouville-Green- Jeffreys- Wentzel-Kramers-Brillouin (L- 
G-J-W-K-B) approximation of the form 

V/ = x V s/n . (172) 
One obtains from the original Dirac equation 

( l7 ^5 + m 4 ) X 4 = (173) 



and 



It follows from (fT75]) that 



7 m V mX 1 = 0. (174) 



detU^d^S + mi) =0. (175) 
Evaluation of the determinant in (I175j) gives the Hamilton-Jacobi equation 

qTtip Sd v S + ml = . (176) 

This shows that the orthogonal trajectories defined by 

m~=g^d^S (177) 

(XT 

are timelike geodesies. The same conclusion follows by applying the L-G-J-W-K-B approximation 
to the second order iterated Dirac equation (|170[) . 
The iterated Dirac equation also gives 

9,SVV = -\{^ 2 S)x l • (178) 

The same result may be obtained by differentiating (|173j) and using (|174| . Thus the spinor is 
parallelly propagated along the timelike geodesies up to direction in spin space. The amplitude 
of the spinor \ l is governed by the expansion of the hypersurface timelike congruence whose 
tangent vector is given by 

vP = —tf w d v S. (179) 

We also have from (|173j) that x % is an eigenspinor of u^^/^. As in flat space it follows that the 
spin tensor 

s = x'tw*' (180) 
x l x l 

satisfies 

S^u v = 0. (181) 

Since x % is parallelly propagated along in direction and since depends only on the direction 
of x % i it follows that the spin tensor is parallelly transported along the timclike congruence, just 
like any other perfect gyroscope. The geodesies are independent of the mass eigenvalue rrii and the 
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polarization state given by \ l - Indeed if the fermion starts off in a given polarization state (with the 
associated mass), it remains in it. In other words, at the L-G-J-W-K-B level, the weak equivalence 
principle, in the form of the universality of free fall, i.e. the statement that all particles fall in the 
same way in a gravitational field independently of their mass, polarization, charge, etc., continues 
to hold. Thus there should be no unusual behavior in the vicinity of a spinning black hole, or indeed 
in the neighborhood of any spinning system due to the Lense-Thirring effect as suggested in [25,26], 
denied in [29] and maintained in [27]. 

6.2 Neutrino mixing matrix 

Here, we briefly review the theory of the neutrino mixing matrix, assuming that the neutrinos are 
Majorana. 

The lepton mixing matrix [30] belongs to the coset U(1) 2 \SU (3), since only phasing of the lepton 
charge eigenstates (f e , v^, v T ), but not the neutrino mass eigenstates V2 1 v^) (which are assumed 
Majorana) is possible. One has 



(182) 
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Thus (y e ,e), (z^,/x), and (v t ,t) are doublets under weak isospin. 
One conventionally fixes the phases so that M takes the form 
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where the three angles #12, $13, and 623 lie in the first quadrant. 

The Jarlskog invariant for the neutrino mixing matrix, defined as in (|10[) but with V now replaced 
by M, is again given by (113p . Note, in particular, that it is independent of the phases a.\ and 0.2- 

Experimentally, parameters of the neutrino mixing matrix are not completely known. According 
to [19], 



sin 2 26» 12 = 0.86; 



-0.03 
-0.04 



0.92 < sin 2 26*23 < 1 



sin 2 20i3 < 0.19, 



(184) 



and there is no experimental information about the Dirac angle 5. Thus, we can certainly deduce 
that there is an upper bound on the Jarlskog invariant for the neutrino mixing matrix, given by 



|J| < 0.049. 



(185) 



For six different neutrino mass eigenstates, as in the seesaw mechanism, a general mixing matrix 
would be an element of U(1) 5 \SU(6), since one would diagonalize a 6 x 6 Hcrmitian matrix. 



6.3 Statistics of J 

We have seen that the parameter space for the neutrino mixing matrix is the six-dimensional single 
quotient U (1) 2 \SU (3) , and that the Jarlskog invariant for the neutrino mixing matrix takes the same 
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form (fTOj) . or ([55]) , as it does for the Kobayashi-Maskawa matrix. Therefore, all results obtained in 
Sees. 13. H and |4. II apply equally to the case of neutrinos. 

For completeness, we quote the results obtained in Sec. 13.11 

(J 2 ) = — w 1.389 x 1CT 3 , (J 4 ) = — - — « 4.960 x 1CT 5 , AJ = —^= w 0.0373 . (186) 
v 7 720 > \ / 201600 ' 12^5 

This can be compared with the experimental bound given in (|185[) . 

One could also repeat the calculations of Sec. 03 assuming particular values for the neutrino 
masses. A strong hierarchy in the neutrino masses would then presumably again lead to "naturally" 
small CP violation from the corresponding mixing matrix. Alternatively, an experimental obser- 
vation of small CP violation for neutrinos would perhaps be an indication of a mass hierarchy in 
neutrinos. At present, neither the magnitude of CP violation nor any values of neutrino masses have 
been measured sufficiently accurately to allow predictions. 



7 Conclusions and Outlook 

In this paper, we analyzed the problem of finding a natural measure on a space of coupling con- 
stants, which in our case was the space of Kobayashi-Maskawa matrices, the double quotient 
U(1) 2 \SU(3)/U(1) 2 . We saw that the measure on this double quotient is nonunique, and we ana- 
lyzed several possible choices of measure on the double quotient. One class of measures was given 
by squashed Kaluza-Klein measures, induced by a Kaluza-Klein reduction of a left-invariant metric 
on the flag manifold. Alternatively, one could take the unique measure on SU(3)/U(1) 2 and simply 
integrate over the left angles. The measure used by Ozsvath and Schucking seemed not to be very 
well motivated from a geometric perspective. 

When calculating expectation values for J, we found that all of the measures we considered led to 
rather similar statistics of J. In each case, the observed value was about three orders of magnitude 
below what one would normally expect; the observed value appears to be finely tuned. The same 
applied to the Ozsvath-Schiicking measure, an extremely squashed Kaluza-Klein measure, or a flat 
measure, which is just the simplest choice and not justified geometrically. 

In Sec. 03 we adopted the different viewpoint that the Kobayashi-Maskawa matrix should not 
be viewed as separate from the quark masses, but that it is really the mass matrices which are 
"chosen" by a yet unknown physical mechanism. We took the observed values for the quark masses 
as an input and chose the simplest distribution which was able to reproduce these observed values, 
while inducing a different measure on the space of Kobayashi-Maskawa matrices. Assuming such a 
distribution, we found that the observed value of J now appears very natural and not finely tuned 
at all. In this statistical approach, regarding the Yukawa couplings determining the mass matrices 
as randomly chosen seems more appropriate than separating quark masses and mixing angles. (On 
submittal of this article to the archive we were informed of an earlier work [31], similar in spirit to 
ours but using different assumptions and methods, which reaches broadly similar conclusions). 
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Our analysis also applies to the case of massive neutrinos, where the predictions will conceivably 
be tested by future experiments. In the standard theory, the Maki-Nakagawa-Sakata matrix [32] 
which appears is naturally an element of the single quotient U(1) 2 \SU(3). Since the right phases 
do not play any role in neutrino oscillations and the relevant J is independent of these phases, the 
calculations are identical to the ones presented here, although with the appropriate values of the /i 
parameters appearing in A and A' . 

In the seesaw mechanism one adds very heavy right-handed neutrinos, and the most general 
mixing matrix would be an element of U(l) \SU(6). This is naturally a Kahler manifold, and the 
measure induced by the Kahler metric can be obtained from the analysis in [17]. We leave a detailed 
treatment of this case, following our approach here, to future work. 

Finally, one could analyze the effects of a fourth generation of quarks on CP violation by repeating 
the calculations for 4 x 4 Hcrmitian matrices. If this generalization spoils the agreement with 
the observed J, one might obtain interesting lower bounds on the masses of a hypothetical extra 
generation of quarks. 
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